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Abstract 

We present the full next-to-leading order (NLO) corrected inclusive cross 
section d^ / dQ"^ / dy / dpT for massive lepton pair production in longitudi- 
nally polarized proton-proton collisions p + p ^ l^l^ Here 'X' denotes 
any inclusive hadronic state and Q represents the invariant mass of the lepton 
pair which has transverse momentum px and rapidity y. All QCD partonic 
subprocesses have been included provided the lepton pair is created by a 
virtual photon, which is a valid approximation for Q < 50 GeV. Like in 
unpolarized proton-proton scattering the dominant subprocess is given by 
Q'(^) + (? — > 7* +' A' so that massive lepton pair production provides us with 
an excellent method to measure the spin density of the gluon. Our calcu- 
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lations are carried out using the method of n-dimensional regularization by 
making a special choice for the 75-matrix. Like in the case of many other 
prescriptions evanescent counter terms appear. They are determined by com- 
puting the NLO coefficient functions for dAa/dQ^ and the polarized cross 
section for Higgs production using both n-dimensional regularization and a 
four dimensional regularization technique in which the 75-matrix is uniquely 
defined. Our calculations reveal that the non-singlet polarized coefficient 
function equals the unpolarized one up to a minus sign. We give predictions 
for double longitudinal spin asymmetry measurements at the RHIC. 

PACS numbers: 12.38.Bx, 12.38.Qk, 13.85.Qk 



2 



1 Introduction 



With the advent of the RHIC at BNL we have a new facihty to study the 
spin structure of the proton, (for a review on the potential of RHIC see [1]), 
which supplements the existing polarized lepton-hadron machines. Polarized 
proton-proton collisions with a very high luminosity and a maximum centre of 
mass energy of ^/s — 500 GeV will provide us with many more details about 
spin distributions than possible with the existing lepton-hadron machines, 
which give very little information about the polarized gluon and sea-quark 
parton densities. An exception is photo-production of heavy flavours where 
the virtual photon is almost on-shell. This process allows us to measure 
the polarized gluon density and has been calculated up to next-to-leading 
order (NLO) in [2]. The gluon density can also be extracted from data on 
prompt photon production, heavy flavour production and massive lepton pair 
production (Drell-Yan or DY process) in proton-proton collisions. The DY 
reaction also provides information on the polarized sea-quark densities. In 
order to measure these densities the above reactions have to be calculated at 
least up to NLO in perturbative QCD. This has been achieved for prompt 
photon production in [3], [4] and for heavy flavour production in [5]. Up to 
the same order calculations have been done for the polarized cross section 
(P /S.a / dQ"^ / dy of the Drell-Yan (DY) process in [6] and [7]. Vector boson 
production in NLO for the cross section di^a jdQ^ has been studied in [8]. 
Notice that these DY cross sections provide us with information about the 
sea-quark density rather than about the gluon density. To determine the lat- 
ter one needs to study the differential distribution d^ /\a / dQ"^ / dy / dpT where 
the transverse momentum px is sufficiently large so that the quark-gluon 
{qg) subprocess dominates the quark- anti- quark {qq) subprocess. According 
to studies in [9] and [10] this already occurs for > (5/2. Notice that the 
quark-gluon subprocess also dominates in prompt photon production at large 
Pt- This reaction has a much larger cross section than the one for the DY 
process, which behaves like at large Q so one might favour it. However 
the former reaction also has its disadvantages, as we have already seen for 
unpolarized prompt photon production, due to experimental and theoretical 
complications with photon isolation criteria. Moreover the DY process has 
an additional large scale represented by the lepton pair mass Q which turns 
out to be a much better scale to use than pt- The NLO QCD corrections 
for the unpolarized DY cross section d^a /dQ^ / dy / dpr were performed in [11] 
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(non-singlet) and in [12] (singlet). A complete NLO QCD calculation for the 
polarized cross section does not exist yet except for the non-singlet part com- 
puted in [13], [14]. One can show, using general arguments, that apart from 
a minus sign the non-singlet DY coefficient function due to the qq-channel 
should be equal to its unpolarized analogue. 

In this paper wc present a complete NLO calculation of the polarized DY 
cross section Aa / dQ"^ / dy / dpT where all partonic subprocesses are included. 
Our findings for the non-singlet coefficient function agree with those obtained 
in [13], [14] and they confirm the relation stated above. In the past most of 
the corrections to physical quantities were performed using the method of 
n-dimcnsional rcgularization [16] to define the singularities which occur in 
loop and phase space integrals. This method was preferred above other reg- 
ularization techniques because it manifestly preserves the Ward-identities in 
non-abelian gauge field theories like QCD. However this advantage gets lost 
when one has to deal with the 75-matrix and the Levi-Civita tensor which 
appear in the clcctroweak standard model and in QCD when one studies 
polarization effects. In this case some Ward identities are violated and have 
to be restored by the addition of so-called evanescent counter terms. There 
exists a large literature about the 75-matrix and the Levi-Civita tensor in 
the context of ultraviolet divergences (sec [15]) but a thorough study on 
how to handle collinear and infrared singularities, which are characteristic of 
partonic cross sections, is still lacking. The most often used prescription to 
compute polarized partonic cross sections is called HVBM after the initials 
of the authors in [16] and [17]. As is shown in [18]- [22] this prescription also 
requires that evanescent counter terms have to be taken into account other- 
wise one gets wrong coefficient functions. In practice it turns out that the 
computation of higher order corrections to physical quantities in the HVBM 
prescription is more complicated then the usual n-dimensional regularization 
technique. This can be attributed to the split up of the n-dimensional space 
in a 4 and an n — 4 dimensional subspace. Accordingly the gamma-matrices 
and the momenta have to split up which complicates the gamma-matrix alge- 
bra and the phase space integrals. To avoid this we introduce in this paper a 
prescription for the 75-matrix which is easier to implement in computer alge- 
bra programs which do not contain a package to deal with gamma-matrices 
following the HVBM prescription. Moreover the n-dimensional phase space 
integrals take their usual form as in calculations where the 75-matrix is not 
present. The evanescent counter terms which we need in our NLO calcula- 
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tion of / dQ'^ / dy / dpT are extracted from an NLO computation of the 
DY polarized differential cross section dAa/dQ'^ and from the Higgs boson 
polarized total cross section. To check our procedure we re-calculate these 
coefficient functions using a four dimensional regularization technique (see 
e.g. [23]) where the 75-matrix is uniquely defined. 

Our paper is organized as follows. In section 2 we introduce our notations 
and discuss extensively the technicalities which are involved if one wants to 
give a prescription for the 75-matrix using n-dimensional regularization. We 
introduce our own definition and determine the corresponding evanescent 
counter terms by comparing our results with those obtained by four dimen- 
sional regularization techniques. In section 3 we present the NLO corrections 
to the differential polarized DY cross section and check the coefficient func- 
tions by re-calculating them in a four dimensional regularization method. 
In section 4 we study the NLO corrections to polarized DY production in 
proton-proton collisions at RHIC and make a comparison with earlier results 
which were obtained in LO only. The long formulae for the soft-plus-virtual 
gluon contributions to the coefficient functions can be found in Appendix A. 



5 



6 



2 Lowest order contributions to the polarized 
Drell-Yan process 

In this paper we consider the semi-inclusive Drell-Yan process 
H,{P,, S,) + H2{P2, S2) ^ 7*(?) +' X' , 

^l+{h) + l-{l2) 

s^{Pi + P2f , = = {h + hf , (2.1) 

where {i — 1,2) represent the incoming polarized hadrons carrying the 
momenta Pj and and spins S^. Further 'X' denotes any inclusive hadronic 
state which is unpolarized. The lepton pair is represented by /"*"/" with 
momenta /i, l2- In this paper we will only consider lepton pairs which have 
a sufficiently small invariant mass Q so that the photon dominates in the 
above reaction and Z-boson exchange effects can be neglected. In addition 
to the kinematical variables in Eq. (2.1) we need the following variables 

T={P,-q)\ U= {P2 - qf , (2.2) 

to obtain the differential cross section 

In the QCD improved parton model the hadronic DY structure function 
dAW^^^'^ is related to the partonic structure function dAWab as follows 

r/^ AW 

xsy^^{s,t,u,Q^^,'). (2.4) 

d t d u 

In the formula above Afa{x, yU^) (a = q, q, g) are the polarized parton proba- 
bility densities where /j, denotes the factorization/renormalization scale and 
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X is the fraction of the hadron momentum carried away by the parton. The 
DY partonic structure function (fAWaj^a2 is computed from the partonic 
subprocess 

ai{Pu Sl) + a2(P2, S2) 7*(?) + ^i(^i) • • • ^m(^m) (2.5) 
and it reads 



m . j3l. / m 



x\AMa,+a,^j*+b,-bj'' , (2.6) 

where i^'aia2 denotes the colour and spin average factor and the polarized 
matrix elements are denoted AM (when we refer to unpolarized structure 
functions, matrix elements and parton densities we drop the A). The partonic 
kincmatical variables are defined in analogy with the hadronic variables in 
Eqs. (2.1), (2.2) as 

s=(pi+P2)', t = {pi-qf, u = {p2-qf. (2.7) 

In the case parton pi emerges from hadron Hi{Pi) and parton p2 emerges 
from hadron H2{P2) we can establish the following relations 

pi = Xi Pi , P2 = X2 P2 , 

s = xiX2S, t = xi{T - Q^) + Q\ u = X2{U - Q^) + Q\ 

- - -^iiT-Q')-Q' 

When parton pi emerges from hadron 7^2(^2) and parton p2 emerges from 
Hi{Pi) one obtains the same expression as in Eq. (2.4) except that T and 
U are interchanged. This result has to be added to Eq. (2.4). When the 

partonic cross section is symmetric under t ^ u one can also use the represen- 
tation in Eq. (2.4) without adding the result where T and U are interchanged 
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provided one makes the replacement A/^^A/^^ _^ A/^^A/^^ _^ A/^^^^/^i. 
Finally note that the relation between the parton probability densities above 
and the parton momentum densities appearing in the parton density sets in 
the literature or PDF libraries, which are denoted by A/^°^(x, //^), is given 
by Af^^^ix,fi')=x AUx,f^'). 

When |AAfaia2p in Eq. (2.6) is calculated up to order «^ one encounters 
four partonic subprocesses which are characterised by the two partons in 
their initial state. In the case of quarks with a mass m ^ they are given 

by 

QiPu si) + q{p2, S2) 7* +' X' , 

I AM,/ = ^ Tr (75/2 ih - m)M-f,^, i^, + m) M^) , (2.9) 
qi{qi){pi,si)+g{p2)^-f* +' X' , 

\AM,,\' = ] e,.xa ^ Tr (m^ 75^ (|5i ± m) M'^t) ^ (2.10) 

?l(?l)(Pl,Sl) +?2(g2)(P2,S2) 7* +' X' 

\AM,,,f = ^ Tr (75^ (^ ± m) M ± m) M^) , (2.11) 

gipi) + g{p2) ^ 7* +' X' , 

1 „A2 70-2 „Al /CTl 

J 6,2.2A2.. ^ e,,^,,,^, ^ Tr (M'^i'^^ M--t) , (2.12) 
4 P2- h Pi - 1' ^ ' 



where M denotes the matrix element which is given by the standard Feynman 
rules. Further the symbol Tr can represent multiple traces when the matrix 
elements are calculated in higher order and for the reaction in Eq. (2.11) one 
can distinguish between — qi and qi 7^ 52- The spin vectors Sj and the 
gauge vectors Zj {i — 1,2) satisfy the properties 



' Pi ) ' 1 ) 1% ' li . 



(2.13) 
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When the (anti-)quark is massless then one has to make the replacements 

75^i {i>i ±m) ^ ±-f5hi^i, (2.14) 

where hi {i = 1, 2) represent the hehcities of the incoming (anti-)quarks and 
the + and — signs on the righthand side hold for the quarks and anti-quarks 
respectively. The definitions above are chosen in such a way that the partonic 
polarized structure function satisfies the property 

^Wa,a, = Wa,a,{+, +) - W^aia.(+, , (2.15) 

with +, — denoting the helicities of the incoming partons. 

The computation of the matrix elements in Eqs. (2.9)-(2.12) and their 
virtual corrections reveals divergences which occur when the momenta over 
which one integrates tend to zero (infrared), infinity (ultraviolet) or collinear 
to another momentum (collinear). The most popular way to regularize these 
singularities is to choose the method of n-dimensional regularization [16] in 
which the space is extended to n dimensions. The singularities are repre- 
sented by pole terms of the type (1/e)^ with n — A + e. This method is very 
useful because it preserves the Ward identities in the case of gauge theories. 
However this is no longer the case when the 75 matrix and the Levi-Civita 
tensor appear like in Eqs. (2.9)-(2.12) or in weak interactions. There is no 
consistent way to generalize these two quantities in n dimensions contrary 
to the ordinary matrix 7^ or the metric tensor Qf^^,. In the literature one 
has proposed various methods to extend the 75 matrix and the Levi-Civita 
tensor to n dimensions but one always needs so-called evanescent counter 
terms to restore the Ward identities. A very popular prescription is the 
HVBM-scheme which was proposed in [16] and generalised in [17]. In this 
approach the gamma-matrices and the momenta have to be split up into a 
4 and an n — 4 dimensional part. Therefore also the integrals over the final 
state momenta have to split up in the same way. Many NLO calculations 
have been done in this scheme (see e.g. [2]-[8], [13], [14]). However this ap- 
proach requires a special procedure to deal with the gamma-matrix algebra 
which is not implemented in the program FORM [24]. Since this program 
is used in our calculations we prefer another prescription for the 75-matrix 
which is given in [25]. It gives the same results as the HVBM-scheme but 
it is much simpler to use in algebraic manipulation programs. Moreover one 
does not have to split up the integrals over the final state momenta and one 
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can simply use the phase space integrals computed for unpolarized reactions 
(sec e.g. [26]). The procedure in [25] is given by 

1. Replace the 75-matrix by 

1 i 

7m 75 = g C/xpar Y Y 01 75 = ^ ^parn Y Y Y Y ■ 

2. Compute all matrix elements in n dimensions. 

3. Evaluate all Feynman integrals and phase space integrals in n-dimensions. 

4. Contract the Levi-Civita tensors in four dimensions after the Feynman 
integrals and phase space integrals are carried out. 

(2.16) 

Note that the contraction in four dimensions only applies to Lorentz indices 
which are present in the Levi-Civita tensors. The last step in (2.16) requires 
that one first has to apply tensorial reduction to all integrals. For simple 
expressions this takes the form 




k^k'' f{k,p,,p2) = Aoo(n)^'''^ + Ai(n)pi'p5: + A22(n)p^p^ 



+ Ai2 (n) p1 P2 + A21 (n) p^ pi , 
j dPS^'^^k^k" f{k,p^,P2) = Boo{n)g^'' + B^,{n)p1p'i + B22{n)p^p'2 

+Bu{n)p1p^,+B2i{n)p^p\, 

(2.17) 

where the coefficients Aij,Bij depend on n = 4 + £. The HVBM-scheme 
or the prescription above automatically reproduces the Adler-Bell-Jackiw 
anomaly (ABJ) [27] but it violates the Ward identity for the non-singlet 
axial vector current and the Adler-Bardeen theorem [28]. In order to obtain 
the correct renormalized quantities one therefore needs to invoke additional 
counter terms [18], which are called evanescent since they do not occur for 
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n = 4. This procedure was used to obtain the NLO anomalous dimensions for 
the spin operators which determine the evolution of the parton spin densities 
see [19] , [20] . Notice that in our case the Drell-Yan process only involves one 
photon exchange. Therefore the 75-matrix does not appear in the virtual 
amplitude so that there is no problem in the ultraviolet sector. Evanescent 
counter terms are also needed in cases where collincar divergences show up 
like in partonic cross sections (see e.g. [22]). This is characteristic for the 
HVBM scheme as well as any other approach similar to the one we propose 
later on. The reason is that in n-dimensional regularization there exists a 
one-to-one correspondence between the ultraviolet divergences occurring in 
partonic operator matrix elements and the collinear divergences appearing in 
partonic cross sections provided both quantities are of twist two type. The 
main problem with each prescription to extend the 75-matrix and the Levi- 
Civita tensor beyond 4 dimensions is that one does not know beforehand 
which Ward identities are violated. Hence in principle all e = n — 4 terms in 
the matrix element might be spurious so that one is never sure to obtain the 
correct result. An example is the spurious terms occuring in the interference 
terms in the quark-quark channel which are discussed below Eq. (3.15). 

A way to avoid this problem is to drop n-dimensional regularization and to 
resort to four dimensional regularization techniques where the 75-matrix anti- 
commutes with the other 7^ matrices and the Levi-Civita tensor is clearly 
defined. However four dimensional regularization methods have their draw- 
backs too because they also violate Ward identities in particular in the ul- 
traviolet sector where it is difficult to apply techniques like Pauli-Villars [29] 
or the introduction of an ultraviolet cut off on the divergent integrals. In 
the case of collinear and infrared divergences four dimensional techniques 
have been applied in a more successful way provided certain rules are re- 
spected. In [23] it was shown that off-shell and on-shell regularization leads 
to the same coefficient functions as those found by using n-dimensional reg- 
ularization as regards the NLO unpolarized DY total cross section. In [30] 
the same result was shown for the process q + q ~^ q + q + 'y* where the 
NNLO coefficient function agrees with the one obtained by the off-shell and 
on-shell techniques apphed in [31]. The off-shell and on-shell regularization 
techniques are defined by (see Eqs. (2.9)-(2.12)) 
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ofF-shell 

Pi ■ Pi = Pi < ^, m = 0, 

Pi represents the quark as well as the gluon momentum , 
on-shell 

Pi ■ Pi — 7^ , Pi is the quark momentum , 

Pi ■ Pi = , Pi is the gluon momentum . (2-18) 

The on-shcll technique has the advantage that the partonic cross sections 
are gauge invariant like the ones computed by n-dimensional regularization, 
where the mass m is set to be zero. The disadvantage is that the collinear and 
infrared divergences due to the gluon have to be regularized by a different 
technique. The off-shell regularization is universal both for gluons and quarks 
but the partonic cross sections depend on the chosen gauge. However the 
same gauge dependence also appears in the operator matrix elements so that 
it cancels in the coefficient functions after subtraction of the former from 
the partonic cross sections. This is checked for the cases mentioned above 
Eq. (2.18). A more serious drawback is that one has to be very cautious 
when collinear and infrared divergences appear together which happens in 
the virtual and soft gluon corrections. As was shown in [23] the off-shell 
assignment as given in Eq. (2.18) leads to the wrong coefficient function 
unless it is computed in the axial gauge (see [32]). Therefore in the case of the 
virtual and soft gluon corrections one has to introduce a more sophisticated 
off-shell regularization which respects the Kinoshita double cutting rules as 
formulated in [33]. In the case when only collinear divergences appear, which 
happens in radiative processes, the off-shell regularization in Eq. (2.18) can 
be successfully applied. 

We will use in parallel two regularization techniques to compute the NLO 
corrections to the double differential structure function in Eq. (2.3). This 
allows us to overcome the above problems associated with all the regular- 
ization techniques and with the evanescent counter terms which have to be 
introduced in the case of n-dimensional regularization. For the first method 
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we adopt a prescription which is shghtly different from HVBM in Eq. (2.16), 
namely 

1. Replace the 75-matrix in by 

i % 
7;. 75 = g e^par Y 7" Y or 15 = ^ eparK Y Y Y Y ■ 

2. Compute all matrix elements in n dimensions. 

3. Evaluate all Feynman integrals and phase space integrals in n dimen- 
sions. 

4. Contract the Lcvi-Civita tensors in four dimensions after the Feynman 
integrals and soft gluon phase space integrals are done. 

5. Contract the Levi-Civita tensors in four dimensions first before the 
hard gluon phase space integrals are carried out 

(2.19) 

Note that by hard gluon phase space integrals we mean all integrals where we 
do not encounter infrared divergences. In the subsequent part of this paper 
we will call this version of n-dimensional rcgularization the RSN-scheme. 
It only differs from the scheme in Eq. (2.16) regarding the procedure in 
contracting the Levi-Civita tensors. We have chosen this approach to avoid 
tensorial reduction of the hard gluon phase space integrals in n dimensions 
as is needed for the scheme in Eq. (2.16) which is equivalent to the HVBM- 
schcmc. For the virtual corrections wc can simply use the method of tensorial 
reduction of the Feynman integrals into scalar integrals which is presented 
in [34] and [35]. Notice that the soft gluon contributions have to be treated 
in the same way as the virtual corrections because they have to be added in 
order to cancel the final state coUinear and infrared divergences in the final 
result. Because there are only a few soft gluon phase space integrals, their 
tensorial reduction is feasible. A drawback of the RSN approach is that one 
has to compute many more evanescent counter terms than in the case of the 
HVBM-scheme. These terms will be extracted from the computation of the 
completely integrated NLO structure function which is one order lower in ag 
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than the NLO correction to the double differential structure function in Eq. 
(2.3). The reason that n- dimensional integration and the contraction of the 
Levi-Civita tensors in four dimensions do not commute can be attributed to 
the metric tensor g^^,^ in Eq. (2.17). When the integration is performed first 
acts as an n-dimensional object which is indicated by the dependence of 
the coefficients of Aij{n) and Bij{n) in Eq. (2.17) on n. When the contraction 
is done first in four dimensions all metric tensors appearing in the matrix 
elements are four dimensional objects provided they have the same Lorentz 
indices in common with the Levi-Civita tensors. This explains the difference 
between the HVBM and the RSN scheme. We checked that in the latter 
scheme |AMqgp in Eq. (2.10) does not depend on the gauge vector I2. This 
independence persists if one adds an unpolarized gluon in the final state with 
gauge vector I to the lower order processes in Eqs. (2.9), (2.10) where the 
sum over the physical polarizations is given by 

P'^^il^p)^ Y: e^(p,«)6^(p,a) = -^'^^ + ^^^4^-^. (2.20) 

a=L,R ' ■ P 

The polarization sum above satisfies the conditions 

p/^*^ = P^'^ = , l^^O. (2.21) 

There is one exception. If there are two polarized gluons in the initial state 
like in reaction (2.12) it turns out that \AMgg\'^ depends on li and I2 so that 
one has to make a specific choice for them. In this paper we will choose 
h = p2 and I2 = Pi- In order to check the RSN method for the computation 
of the 2 — 3 body hard gluon processes we compute the double differential 
partonic structure function in four dimensions using the off-shell regulariza- 
tion technique in Eq. (2.18). This will be our second regularization method. 

In the final part of this section we have to compute the evanescent counter 
terms which show up in the RSN-scheme. They arc extracted from the NLO 
totally integrated structure function in Eq. (2.6) which has to be computed 
using several regularization schemes. Here we assume that only four dimen- 
sional regularization techniques provide us with the proper answer since the 
75-matrix anti-commutes with the other 7^-matrices. On the Born level we 
have the following partonic subprocesses (see Eq. (2.5)) 



q{pi) + q{P2) ^ 7* , 9{Pi) + 9{P2) ^ H , (2.22) 
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where H represents a scalar particle e.g. the Higgs boson which like the 
photon has a mass indicated by the same variable Q. The reason that we 
need the latter process is that the totally integrated DY structure function 
in NLO only provides us with information about the splitting functions AP^^ 
and APqg but we also need to know AP^g and AP^^ which can be extracted 
from scalar boson production. In lowest order the matrix elements are 



I am: 



(0) 



-Ns 



1-'- 
2 



I AM. 



(0) 



99- 



--{N'-l)s\ (2.23) 



where denotes the number of colours. The reason that the latter matrix 
element has one power in s more than the former can be attributed to the 
effective scalar-gluon-gluon coupling which has a dimension which is equal 
to the inverse of a mass. This coupling like the electromagnetic coupling in 



I AM 



(0) 



qq^J* I 



has not been included in the definition of AW. On the Born 



level the polarized structure functions become 



AW, 



(0) 



QQ 



99 



1 



m - 1 



5(1 — x) , X 



(2.24) 



Notice that for the DY process we have divided by 1 — ^/2 whereas for scalar 
particle production we have removed an overall term Q^/S in order to get 
unity on the right hand side of AWj^ and AM/j^-' respectively. In NLO the 
matrix elements in the RSN scheme as well as in the off-shell mass scheme 
read 

q{pi) + q{p2) g{ki) + 7*(g) , 



lAM^y J 2 



I at: 



(1)|2 



QQ 



-NCpg^lAT^^ , 
I tu 



— e 




2 Q I 2 w 

Pi 75" +P2 — 



Kqq = 



1 



(2.25) 
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(7V2-1)T^/|AtS)|2 



4^ + 2- -2- 

t s s t 



19 



N {N^ - 1) ' 

q{q){Pi) + 9{P2) q{q){ki) + H{q) , 



(2.26) 



|AM, 



gq^q H\ 



1 



NCpg' 



u 



2Q' 
Pi — 



K. 



9<1 



|AM, 



gg^g H\ 



N {N^ - 1) ' 

9{Pi) + 9{P2) ^ g{ki) + H{q) , 



(2.27) 



^N{N'-l)g' 



2s(s- g2) +^ ^ 2s(s - + 4 



+— + ^ — + 

tu s s 



+Pi — +P2 



K. 



99 



(2.28) 



Here g represents the strong coupling constant and the colour factors are 
given by 



2N 



(2.29) 
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Notice that one has to put = in the equations above if the matrix 
elements are understood to be used in the RSN scheme where the quarks 
are massless and the phase space integrals are evaluated in n dimensions. In 
the case of the off-shell regularization method one has to put s = 0. Also 
we have neglected higher powers in p| since they do not contribute to the 
structure function in the limit pf 0. Further the off-shcll matrix elements 
are gauge dependent and we have taken l'^ = = for the polarization sum 
in Eq. (2.19) (Feynman gauge). In the case of the on-shell mass assignment 
we obtain 



I AM 



(1) 



{t — m?) {u — m?) 



+ 4:m 



s{t — rin?Y 



s {u — w?Y 



(2.30) 



I AM 



(1) 



99^9 7 



(Ar2-i)T// 



t — m? s — m? 



I AM 



(1) 



gq^q H\ 



—m 



4g^ 



(t — m?y 
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(2.32) 



Notice that we have omitted the on-shell matrix element lAM"^^^-*!^ because 
the gluon is massless so there is no regulator mass for the coUinear divergences 
when we integrate over the final state momenta. 

Let us first start with the computation of the structure functions in the 
off-shell and on-shell mass assignments because they immediately lead to the 
correct coefficient functions in the MS-scheme. Integration over the final state 
momenta k-i and q provides us with the following results. For the off-shell 
mass assignment we get 
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x)'^ - 12 + 16x 



(2.36) 



where we have introduced the shorthand notation 

,2 



etc 



47r (47r)^ 



(2.37) 



Notice that we have omitted the soft-plus-virtual gluon corrections in Eq. 
(2.33) and Eq. (2.36) which are proportional to 5(1 — x). As was pointed out 
in [23] these terms cannot be correctly computed using the off-shell assign- 
ment according to Eq. (2.18). In the case of the on-shell mass assignment 
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we can compute the latter correction and we obtain 
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(2.39) 
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(2.40) 



Eqs. (2.34) and (2.39) arc in agreement with [36]. A comparison with the 



unpolarized structure functions Wg^^ and Wg^' as computed for instance in 
[23] reveals the following relations 
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The last term is due to chiral symmetry breaking since the quark has a mass 
m 7^ 0. 

We will now present the coefficient functions following from mass factor- 
ization denoted by 



k,l=q,g 



(2.42) 
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In the formula above all collinear divergences are absorbed by the kernels 
ATki which depend on In/x^/p^, In/x^/m^ or 1/e depending on the regular- 
ization method used. Here /i represents the factorization scale which also 
enters the finite DY coefficient functions AWki. In NLO the latter are given 

by 
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qq 
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99 
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(2.43) 

In the case of the off- and on-shell regularization the kernels are given by the 
operator matrix elements as follows 

Ar,, = AA,, , AFg, = i AAgg , AF,, = AA^^ , AF,, = AA,, . (2.44) 

For the off-shell mass assignment the latter can be found in [19], [21] and 
they are given by 

AA. = AAf + AAf , 

AA. = AAf«^^ i = (2.45) 

The term A^^*-**^ is characteristic of the off-shell mass assignment and it 
vanishes when the equations of motions (EOM) are applied to the external 
quark indicated by the momentum p so that it does not show up in the on- 
shell case. In the Feynman gauge and renormalized in the MS-scheme they 
read 
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(2.47) 
, (2.48) 
, (2.49) 
(2.50) 
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The results above are calculated using the HVBM prescription according to 
Eq. (2.16). In this case the renormalization constants become equal to 



7HVBM 
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where AP-j^ are the lowest order polarized sphtting functions (see e.g. [19], 
[20]) given by 

AP(J) = r/[16x-8], 
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(2.53) 



However there is one exception in the case of Zqq. For the HVBM prescription 
the latter is given by 
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The last term in is the evanescent counter term which is necessary to 

ensure that the non-singlet axial vector current does not get renormalized. 
This implies that the first moment of ^A^^^^ must be unity in all orders 
of perturbation theory. Notice that in higher order the evanescent counter 
term even becomes infinite (see [18], [21]). ^ Using the same counter terms as 
in Eqs. (2.52), (2.54) one obtains the following expressions for the operator 
matrix elements in the on-shell mass assignment. Using the notation in Eq. 
(2.45) the operator matrix elements /^A^J^^^ become equal to 



5{l-x) + asCF 



l-x 



- 2 - 2x I In At 
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-2 + 6x + (5(1 - x) i 3 In At + 4 



(2.55) 



^It appears that in the case of the naive 75-matrix prescription, where this matrix 
anti-conimutcs with the other 7-matrices in n dimensions, an evanescent counter term is 
not necessary in the non-singlet case. 
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qg 



8a; - 4 > In 



(2.56) 
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(2.57) 



and A^f = 0. Notice that the first moment of AA°^ in Eq. (2.55) is still 
not equal to unity even after subtraction of the evanescent counter term. This 
is because the non-singlet axial vector current is not conserved for massive 
quarks and therefore has to undergo a finite renormalization. This additional 
term, which is equal to ttgCp A{1 — x), will compensate the last term in Eq. 
(2.41) while computing the difference in Eq. (2.43). After mass factorization 
one obtains the same coefficient functions irrespective whether the off-shell 
or the on-shell regularization is used. The results are 
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Notice that the coefficient of tlie b(\ — x) term in Eq. (2.58) can only be 
obtained from tlie on-sliell sclieme. Tlie latter scheme could not be applied 
to the process gf + gi— >gi + i7in Eq. (2.28) so that the 5(1 — a;)-term in Eq. 
(2.61) can only be inferred from n-dimensional regularization as HVBM or 
RSN. Furthermore from Eq. (2.58) one derives 
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(2.62) 



This relation, which only holds for coefficient functions, has to be satisfied 
irrespective which regularization scheme is chosen. 

Now we want to compute the coefficients above using the RSN approach 
presented in Eq. (2.19) and determine the corresponding evanescent counter 
terms. The DY partonic structure functions are given by 
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This scheme violates the relation in Eq. (2.62) in a similar way as the on-shell 
result in Eq. (2.41). This is because the n-dimensional extension of the 75- 
matrix breaks the chiral symmetry. However this is unphysical because the 
quarks are still massless. The kernels for the RSN prescription follow from the 
requirement that the RSN-scheme must lead to the same coefficient functions 
as the four dimensional rcgularization methods. From mass factorization (see 
Eq. (2.43)) and the coefficient functions presented in Eqs. (2.58)-(2.61) we 
infer the kernels 
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where the additional terms are characteristic of the RSN-scheme. 
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Finally we want to show that the prescription given in Eq. (2.16) leads 
to the same results as the conventional HVBM approach in [16] and [17]. A 
straightforward calculation gives 
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(2.72) 

(2.73) 
(2.74) 



The results in Eqs. (2.71) and (2.72) agree with those in [6] which were 
obtained by the conventional HVBM approach. The kernels for the HVBM 
prescription are obtained according to the same requirement as used for the 
RSN-approach above Eq. (2.67) and they read 
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As in the case of unpolarized quantities the kernels above are equal to 
the operator renormalization constants in Eqs. (2.52), (2.53), except that 
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^pHVBM AZ^v^^ (sec also Eq. (2.44)), provided the same 75-matrix 

prescription is used for ultraviolet and collinear divergences. Further notice 
that RSN and HVBM schemes only affect the regular part of the kernels and 
do not alter the singular part represented by 5(1 — x). 

Summarizing the results obtained in this section we have introduced a 
new scheme, called RSN, which is designed to regularize the hard gluon con- 
tributions to the double differential structure functions in Eq. (2.4) when the 
latter are computed in NLO. For that purpose we calculated the evanescent 
counter terms by the requirement that in the MS scheme one should obtain 
the same coefficient functions as obtained by four dimensional regularization 
methods. One of the most important results is the relation between the po- 
larized and unpolarized structure functions in Eq. (2.62) which has to be 
satisfied irrespective of the prescription used for the 75-matrix. That this 
relation has to emerge from the calculations can be simply inferred from the 
off-shell regularization technique (see the result for Wqq in Eq. (2.41)). In 
four dimensions the two 75-matrices in the non-singlet part of the matrix 
element in Eq. (2.9) can be brought together since they are present in the 
same trace so that they yield unity. In this way the matrix element becomes 
the same as for unpolarized reactions. 
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3 NLO corrections to polarized lepton-pair 
production 

After having specified the RSN-prescription for the 75 matrix and the Levi- 
Civita tensor in the previous section we present the NLO corrections to the 
quantity 
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Here r{x) denotes the gamma-function and the spherical factor is defined 

by 



'5£ = expQ(7E-ln47r)^ . 



(3.2) 



Notice that the expression in Eq. (3.1) holds for the Born reactions in Eq. 
(2.25), (2.26) as well as for the virtual corrections presented below. The 
calculation proceeds in the same way as e.g. done in the case of unpolarized 
processes hke heavy flavour production [37] or Higgs production [26] . Starting 
with the virtual contributions we compute the one-loop corrections to the 
reactions in Eqs. (2.25), (2.26). The Feynman graphs can be found in Fig 5 
of [30]. First we perform tensorial reduction of the Feynman integrals into 
scalar integrals following the procedure in [34] , [35] . The scalar integrals can 
be found in Appendix D of [38] and we use n-dimensional regularization for 
the ultraviolet, infrared and coUinear divergences. Finally we contract the 
Levi-Civita tensors in four dimensions. The results which hold in the HVBM- 
as well as RSN-schemes are given by 
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Apart from an overall minus sign the above result agrees with the unpolarized 
virtual contribution in [11]. Note that e terms in lAT^^"*!^ (see Eq. (2.25)) 

differ from those in the unpolarized case, denoted by which causes 

the breakdown of the relation in Eq. (2.62). Also we find 
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Here the function Li2(a;) in the expressions above denotes the di- logarithm 
which e.g. can be found in [39]. Next we compute all 2 — > 3-body processes 
which contribute to the structure function 

^ AW^ll ^ J_ ^ SI (tu-Q^ s,y/' ( s, \ 
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with S4 = s + t + u — Q , (3.5) 



32 



where | AM^^),^^^^^^ ^, ^ is the second order matrix element integrated over 
the polar angle 9i and the azimuthal angle 62 so 



— >6i62 7* I 

= /^'^ d^i (sin 6,)'+^ £ de^ (sin ^2)^ \AM^%,_^,^b. r (^i' ^2) 1' 

^ / d^ln-l|AMi^)^_,,^,^ ^.(^1, ^2)r ■ (3.6) 

The Feynman graphs for the 2 — > 3 body parton subprocesses can be found 
in Figs. 6-9 in [30] . These reactions have been calculated using two different 
regularization methods for the collinear divergences. In n dimensions we use 
the RSN prescription mentioned in Eq. (2.19). In this case the collinear 
divergences in Eq. (3.6) manifest themselves as pole terms of the type 1/e. 
For n = 4 we adopted the off shell regularization method in Eq. (2.18) where 
the collinear divergences show up as logarithms Ins/p^. Notice that the off- 
shell regularization technique is gauge dependent and we have computed Eq. 
(3.6) in the Feynman gauge (see the comments below Eq. (2.29)). The gauge 
dependent terms will be cancelled after mass factorization via the operator 
matrix elements in Eqs. (2.45)-(2.50) which are computed in the same gauge. 
Decomposing the matrix element in Eq. (3.6) into colour factors one obtains 
the following results for the non-singlet qq and the singlet qg subprocesses 
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(3.8) 



The colour and spin average factors Kqq and Kqg are given in Eqs. (2.25) 
and (2.26) respectively. Further the indices Cp and Ca refer to the highest 
power of the colour factor multiplying the corresponding matrix elements. 
Finally we have also explicitly indicated the charge of the quark indicated 
by Cq. For the lowest order matrix elements in section 2 we have implicitly 
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taken it into account. However as we will see below in higher order processes 
the Feynman graphs contain different species of quarks. Their charges will 
appear in the corresponding coefficient functions and they indicate how they 
have to be combined with the parton densities. In NLO we encounter some 
new subprocesses. The first one is given by quark-quark scattering (non- 
identical and identical quarks) represented by 

gi + 52 ^ gi + g2 + 7* , ?i 7^ ^2 , 
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In Eqs. (3.9), (3.10) C, E and F refer to the graphs in Figs. 8, 9 in [30]. 
The second subprocess is quark-anti-quark scattering 
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In Eqs. (3.11), (3.12) and (3.13) A, B, C and D refer to the graphs in Figs. 
7, 8 in [30] . The last new process which shows up in NLO is given by 
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with the same notation as explained below Eq. (3.8). If the matrix elements 
are computed in n dimensions one has to be careful with the interference 
terms indicated by the indices CD, CDEF in Eq. (3.10) and ABCD in 
Eq. (3.13). These terms do not lead to collinear divergences so they can 
be computed in four dimensions. However if they are computed in n dimen- 
sions the prescription for the 75-matrix leads to additional terms which are 
proportional e = n — 4. After integration over the angles in Eq. (3.6) these 
terms contribute because they are multiplied by pole terms 1/e. Since a finite 
quantity cannot depend on the regularization these terms are considered to 
be spurious and they should be discarded. Note that without the 75-matrix 
the terms proportional to e cancel automatically among themselves. Some 
of the expressions i.e. Eqs. (3.9), (3.10), (3.12), (3.13) are singular in the 
limit S4 ^ 0, where S4 is defined in Eq. (3.5), which is due to soft gluon ra- 
diation or soft (collinear) fermion pair production. This will lead to infrared 
singularities when the differential cross section in Eq. (3.5) is convoluted 
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with the parton densities. Therefore the cross section will be split up into a 
hard gluon (54 > A) part and a soft gluon (54 < A) part (see e.g. [26]). The 
latter is defined by 



Only the singular part of | AM^^^^^jj^j,^ 7* which behaves as I/S4, con- 
tributes to the above integral whereas the non-singular terms vanish in the 
hmit S4 — > 0. The calculation of the soft gluon cross section has to proceed in 
the same way as done for the virtual corrections in Eqs. (3.3), (3.4). First we 
apply the eikonal approximation to the singular part of the matrix elements 
and apply tensorial reduction to the soft gluon phase space integrals before 
we contract the Levi-Civita tensors in four dimensions. The results which 
hold in the HVBM- as well as RSN-schemes are given by 
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xIAT, 
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We have added the soft ghion contributions from Eqs. (3.9), (3.12) and 
(3.13) into the above expression. Like in the case of the virtual contribution 
in Eq. (3.3) the expression above is, apart from an overall minus sign and the 
£-dependence of lAT^^^-*!^, equal to the unpolarized structure function. Note 
that this £-dependence leads to the breakdown of the relation in Eq. (2.62). 
Prom the reaction in Eq. (3.9) we obtain 
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|at(i)|2 
\ qg I 



(3.18) 



The infrared and final state collinear divergences are cancelled upon adding 
the virtual contributions in Eqs. (3.3) and (3.4) to the soft gluon contribu- 
tions above. The ultraviolet divergences are removed by subtracting from 
d'^ AW^fl^/dt du (ai, 02 = q, g) the counter term 



(1) 



2/?o d'^Wl, 



^ £ dtdu 



a2 



With f3o = ^CA-^Tfnf, (3.19) 



where we have chosen the MS renormalization scheme. Here denotes 
the lowest order contribution to the /3-function. After renormalization one 
still has to perform mass factorization to remove the remaining collinear 
divergences. This is achieved by the formula 



d^ AW 



a\a2 



dtdu 



{s,t,u,e) = I dxi_ I ^^r^^^^(2.^^£-)^r^_^^^(2.2,£) 



bi,b2=q,9 
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s, t, u) 



bi,b2=Q,9 



61 ai 



6202 ® s 



d^ ^W^,b, 
dtdu 



, (3.20) 



with the following definitions 



s = xiX2S, i ^ xi{t - Q"^) + 



u 



X2{u-Q^) + Q\ (3.21) 



In Eq. (3.20) d'^ A.Waia2 represents the singular structure functions which 
contain the collincar divergences indicated by e. These divergences are re- 
moved by the kernels AF^^a^ leaving the finite coefficient functions d"^ AWbib2- 
For the RSN prescription the kernels are given in Eqs. (2.67)- (2.70). In the 
case of the off-shell mass assignment they are represented by the operator 
matrix elements in Eqs. (2.45)- (2.51). Both the kernels and the coefficient 
functions depend on the mass factorization scale /j,. Both regularization tech- 
niques lead to the same coefficient functions renormalized in MS-scheme. The 
coefficient functions originating from the soft-plus- virtual gluon contributions 
are sufficiently short that they can be published and one can find them in 
Appendix A. The hard gluon parts are too long to be published. They exist 
as FORM [24] files and they are available on request. Here we can only show 
those parts which behave like I/S4. In the soft gluon limit where S4 ^ they 
read 



d"^ AiyK>^° 
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s4-*o dtdu 
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|AMW|2. (3.23) 

Another important result which emerges from our calculation is that we find 
the same relation as in Eq. (2.62) for the non-singlet part of the coefficient 
function i.e. 

at a u at a u 

which holds for the hard and soft-plus-virtual part independently (see also 
the comment below Eq. (A.l). The sum of these two parts is in agreement 
with Eq. (60) in [14] provided one omits the interference contribution called 
2(74 which originates from the non-singlet part of the gg-channel. The relation 
in Eq. (3.24) follows from chiral symmetry because the quarks are massless. 
Therefore the finite coefficient functions should respect this relation irrespec- 
tive which regularization method is used. Finally we want to comment on 
the scale /x. In the computation of the radiative corrections we have assumed 
that the renormalization scale //^ is equal to the mass factorization scale 
If one wants to distinguish between these scales one has to substitute 

a^iii^) = a.iixl) [l + as{iil)(3o In ^\ , (3.25) 



in all finite expressions. 
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4 Differential distributions for the process 

p + p^7* V X' 

In this section we present the differential distributions in Eqs. (2.3), (2.4) 
for lepton-pair production in proton-proton coUisions at the RHIC and make 
a comparison with similar results in previous work. In practice one is not 
interested in the distributions which depend on T and U in Eq. (2.2) but in 
the differential cross section given by 

^3 A^PP ^3 A^PP 

iSy^^y^Q') = ASQpT Z , rr iS,T,U,Q') , (4.1) 



dQ dpT dy^ '"^''"'^ ' ^"^'dQ^dTdU 

where y and pt denote the rapidity and transverse momentum respectively. 
Neglecting the masses of the incoming hadrons we have the following relations 



T = - ^/S ^pj, + Q^ coshy + v^Pt + sinhl/, 

U = Q^ -VS^p^ + Q^ coshy-yfS^pl + sinh|/ . (4.2) 
The kinematical boundaries are 

g2_5<T<0, -S-T + Q^ <U <^^-^ + Q\ (4.3) 

from which one can derive 

with p-^ ^i-^^'X (4-4) 



4 5* cosh y 



or 



(S - Q2)2 

-I/max <y< I/max , < < — = PT,rn^ , 



1 l + yT^ 4 ^ (p^ + Q^) 

With ym^^-ln ^_^^_^^ , ^^--O^TQ^- ^'-'^ 
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Since the cross section in Eq. (4.1) diverges for Pt ^ we cannot perform 
the integral over this kinematical variable down to zero. However the full 
integration over the rapidity can be carried out and one obtains the transverse 
momentum distribution 

d Q dpT J-yma^ d Q d pt d y 

with i/max given in Eq. (4.5). Finally we define what we mean by leading order 
(LO) and next-to- leading order (NLO). In LO the differential cross section 
is determined by the leading logarithmic approximations to the couphng 
constant asi^J'r) in Eq- (3.25) and the polarized parton densities Af^^{x. /i^) 
in Eq. (2.4). In NLO the latter two quantities are replaced by their next-to- 
leading logarithmic approximations. 

In the subsequent part of this section we will study the dependence of 
the cross sections as defined in Eqs. (4.1), (4.6) on input parameters like the 
QCD scale A, the factorization scale n and the dependence on the chosen set 
of polarized parton densities. Notice that we have adopted = ^ for the 
renormalization scale. In our computations the number of light flavours is 
taken to be n/ = 4 which holds for the running coupling, the DY coefficient 
functions and the polarized parton densities. Further we have chosen for 
our plots the polarized parton densities provided by GRSV [40] (here called 
GRSVOl) and BB [41], which were determined in the MS renormalization 
scheme. Notice that these sets do not contain charm and bottom quark 
densities. The densities of GRSV standard scenario were constructed to fit 
the available data together with positivity requirements and to satisfy two 
SU(3) flavour group sum riilcs. The GRSV collaboration have also given 
another set of parton densities, the so-called valence set, where the sum 
rules hold for the SU(2) flavour group as suggested originally in [42]. It 
turns out that in the standard scenario the gluon density is larger than in 
the valence scenario so that the DY cross section will be dominated by the 
qg process in the former scenario provided the transverse momentum will 
be sufficiently large. However in the valence scenario the sea-quark density 
becomes important too. Therefore there will be a competition between the 
qg ^ and qq subprocesses and it will be very hard to disentangle these two 

^Note that the notation qg represents the subprocesses qg and qg. The same holds for 
the notation qq representing the reactions qq and qq. 
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densities. The BB collaboration have presented two analyses of polarized 
deep inelastic scattering data. At small x the densities behave as x"" where 
the ac for the gluon and the as for the sea quark are related by = og + c. 
In the first scenario BB assume c = 0.9 in LO (ISET=1) and c = 1.0 in NLO 
(ISET=3). In the second scenario BB assume c = 0.6 in LO (ISET=2) and 
c = 0.5 in NLO (ISET=4). All sets above arc presented in LO and NLO 
with the A4 and the corresponding values for as{Mz) in Table 1. Besides the 
polarized cross sections dAcr we also want to compute the double longitudinal 
spin asymmetry defined by 

Al.^^, (4.7) 

where da denotes the unpolarized cross section which is calculated up to NLO 
in [12]. ^ In order to compute the latter we adopt the GRV [43] parton density 
set (here called GRV98) which docs not contain any charm or bottom quark 
densities cither (for more details sec Table. 1). Before presenting the results 
we note that our computer programs were checked by reproducing Fig.4 in 
[14] for the LO polarized and unpolarized nonsinglet pT-distributions. The 
difference between their results and ours is numerically small which might 
be due to a different choice of a^. Notice that the interference term 2(j4 
mentioned below Eq. (3.24) is completely negligible. Also we reproduced 
several results in [10] including in particular Fig.3 for the unpolarized NLO 
inclusive pr distribution and Fig. 8(b) for the LO contribution to v4ll(pt) 
in Eq. (4.7). 

For our plots we choose y/S — 200 GeV which is a representative CM. 
energy for proton-proton collisions at the RHIC. The factorization scale is set 
/Li^ = +Pt unless mentioned otherwise. We begin with Figs, la, lb which 
show the pt dependence in LO and NLO of Eq. (4.6) for Q = 6 GeV/c with 
the set GRSVOl (standard scenario). Both in LO and NLO the qg subprocess 
is positive but the qq channel yields larger contributions than the qg channel 
when Pt < 3 GcV/c. For larger transverse momenta the qg subprocess 
dominates and the qq contribution drops off rapidly where it even becomes 
negative whenpT > 27 GeV/c. At ~ 6 GeV/c the cross section due to the 
qq subprocess is only one quarter of the contribution of the qg channel. Hence 

^We re-calculated the unpolarized DY cross section and found agreement with [12]. 
However we used our versions of the coefficient functions for the plots in this paper. 
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GRSVOl (LO, standard scenario) 




175 MeV 


«r(Mz) 


= 0.121 


GRSVOl (NLO, standard scenario) 




257 MeV 


af^O(Mz) 


= 0.109 


GRSVOl (LO, valence scenario) 


A^° = 


175 MeV 




= 0.121 


GRSVOl (NLO, valence scenario) 




257 MeV 


«f^«(M^) 


= 0.109 


BB (LO, scenario 1) 


ALO = 


203 MeV 


a^O(M^) 


= 0.123 


BB (NLO, scenario 1) 




235 MeV 


«r^(M^) 


= 0.107 


BB (LO, scenario 2) 


A^ = 


195 MeV 


«f^O(Mz) 


= 0.123 


BB (NLO, scenario 2) 




240 MeV 


af^°(M^) 


= 0.107 


GRV98 (LO) 


A^o = 


175 MeV 


«f^°(Mz) 


= 0.121 


GRV98(NLO) 


A^LO = 


257 MeV 


ar°(Mz) 


= 0.109 



Table 1: Polarized and unpolarized parton density sets with the values for 
the QCD scale A4 and the running coupling as{Mz). 

we see that when pt > Q/2 the qg channel begins to dominate over the other 
subprocesses. The NLO results from all subprocesses are presented in Fig. lb 
where we note that the contributions from the gg and qq channels are small 
and negative so that we plot their absolute values. The possibility that the 
polarized gluon density can be measured in the DY process at large pt has 
already been stressed by several authors using LO perturbation expressions, 
for recent work see [10]. We have now demonstrated that this conclusion is 
unaltered when the NLO contributions are included. 

In Figs. 2a, 2b we do the same as in Figs, la ,1b but now the plots 
are presented for the GRSVOl valence scenario. In this case the qq channel 
also dominates the cross section at small transverse momenta but this con- 
tribution is negative over the whole pt range so we plot its absolute value. 
On the other hand the qg subprocess only yields positive contributions. The 
net effect is that below ~ 5 GeV/c the LO polarized cross section due 
to both contributions is negative because the qq reaction is dominant in this 
region. Forpr > 5 GeV/c the LO cross section becomes positive since in this 
region the ^gi-channel is more important. Hence we plot the absolute value 
of LO(sum) in Fig. 2a. This overall picture is not altered by including the 
NLO corrections so also the NLO (sum) is negative for small px and positive 
for large pt and we add the plot of its absolute value to Fig. 2a. In Fig. 2b we 
show the individual NLO contributions. The gg and qq channels are small 
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and negative so we have plotted their absolute values. Therefore the main 
difference between the standard and valence scenarios shows up in the region 
Pt < 5 GeV/c where the standard scenario provides us with a positive NLO 
distribution while the valence scenario yields a negative one. Note that there 
is nothing unphysical about this result because the polarized cross section 
is the difference between two helicity projections. If the valence scenario is 
correct one will need high statistics in the region where the Pt distribution 
changes sign, say for 4 < < 6 GeV/c. 

Given the changes in signs in the above results it is clear that rapid- 
ity plots integrated over a range in p^ will depend strongly on the range 
available experimentally and will also vary dramatically . Therefore we will 
only present the rapidity distributions in Eq. (4.1), at a fixed representative 
Pt- In Figs. 3a, 3b we show the LO and NLO corrected distributions with 
Q — 6 GeV/c and pt = 2 GeV/c using the GRSVOl standard scenario par- 
ton densities. In this region the qq channel dominates the cross section and 
it leads to a positive contribution in the whole rapidity range. Since the qg 
channel has two maxima at symmetric values of y the total LO distribution 
has a slight minimum at y = which is enhanced when the NLO corrections 
are taken into account. This is explained by Fig. 3b where we see that the 
qg, gg and qq contributions are all negative near y = so the dip is more pro- 
nounced. At larger pr the relative contributions from these channels changes 
but the qg becomes a larger fraction of the total so that the double peaked 
nature of the NLO distribution will remain. We repeat these two plots for 
the valence scenario parton densities in Figs. 4a and 4b. However we have 
now chosen p^ = Q GeV/c where the both the LO and NLO cross sections 
are positive. Here one clearly sees that the rapidity distributions are nega- 
tive near y = 0, where they have pronounced minima. This is wholly due 
to the qq subprocess which yields negative contributions to the cross section 
over the whole rapidity range. The qg subprocess is positive for all values 
of the rapidity and compensates the qq contribution leaving a much smaller 
cross section than shown by the standard scenario in Figs. 3a, 3b. This 
was already expected from the transverse momentum distributions in Figs. 
2a,2b where the cross section decreases at small values of pr and is actually 
negative for px < 5 GeV/c. At larger px the double peaked nature of this 
distribution will prevail as the relative contribution from the qg channel in- 
creases. Note that the integrals over the above rapidity distributions agree 
with the values of the pr distributions in Figs, la, lb and Figs. 2a, 2b at 



44 



Pt = 2 and px = 6 GeV/c respectively. 

In Fig. 5 we have shown the LO and NLO corrected cross sections for the 
GRSVOl standard scenario at different values for the di-lepton pair invariant 
mass Q. When we increase the value of Q the pt distributions decrease in 
LO as well as in NLO. Further we observe that the NLO corrections become 
larger when the transverse momentum increases. This will become more 
clear when we show the A'-factor later on. The dependence on Q is also 
investigated for the ^-distributions shown in Fig. 6. Here the the NLO 
correction is conspicuous at Q = 2 GeV/c whereas at larger values of Q 
there are hardly any differences between the LO and NLO corrected cross 
sections. 

Radiative corrections are marred by theoretical uncertainties. The first 
one concerns the choice of the factorization and renormalization scale 
Since the perturbation series is truncated at a certain order of ag physical 
quantities depend on yU. Notice that this dependence will disappear when all 
orders are taken into account. The sensitivity to the scale above is exhibited 
by the ratio 



If this ratio is close to unity and almost independent of /i the higher order 
corrected cross section will be very reliable. In Fig. 7a we use the GRSVOl 
standard scenario and investigate the ratio in Eq. (4.8) for the cross section 
in Eq. (4.6). Here we choose as central value — Q'^ + p^ and we vary fi 
from 0.2 yUo to 5/io at Q = 8 GcV/c. Note that the scale on the fi/fio-axis is 
logarithmic. For the transverse momenta we take Pt = 2, 10 and 20 GeV/c. 
In Fig. 7a the LO corrected N{pt, Ai/z^o) are the upper curves for small nf/J^o 
whereas the lower curves are the NLO corrected N{pt, n/fJ-o)- One sees that 
the scale dependence diminishes in going from LO to NLO which indicates 
better predictive power in NLO perturbation theory. The analogous plot is 
shown in Fig. 7b for d^Aa{pt,y,fi)/dQdpxdy as a function of y, with y = 
and y = ±1 where now IJ^q = + Pt and pt = 2 GeV/c. Since the |/-plots 
are symmetric the values at y = 1 are identical to those for y = —1 so there 
are only four curves in Fig. 7b. Again the LO curves are above the NLO ones 
for small ji/jJiQ. Like in the case of the transverse momentum distributions 
there is an improvement in the scale dependence in going from LO to NLO 
perturbation theory. 




(4.8) 
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Besides the dependence on the factorization/renormahzation scales there 
are two other uncertainties which affect the predictive power of the theo- 
retical cross section. The first one concerns the rate of convergence of the 
perturbation series which is indicated by the X-factor defined by 

This quantity is plotted as a function of pr in Fig. 8a for (5 = 2,6 and 
10 GeV/c in the GRSV standard scenario and we see that it varies from 
approximately 1.2 at = 2 GeV/c to 1.7 at pt = 30 GeV/c. The corre- 
sponding rapidity plots for fixed Pt = 2 GeV/ c in Fig. 8b vary from 0.8 to 1.5 
for the case of Q = 2 GeV/c and from 1.2 to 1.5 for the case oi Q = 6 GeV/c. 
Notice the dips around y — 0. They are explained in the discussion of Figs. 
3a, 3b. The minimum in the LO rapidity distribution becomes deeper when 
the NLO corrections, mainly coming from the qg subprocess, are included. 
The study of the iC-factors reveal that the NLO corrections are appreciable 
(sometimes more than 50%) except at small pt where fortunately the cross 
section attains its maximum. 

A third uncertainty which has to be solved by experiment are the parton 
densities originating from different parametrizations. To study the depen- 
dence on these densities it is better to reduce the effect of the choice of 
factorization scale and the various i^-factors. This is achieved by plotting 
the double longitudinal spin asymmetry as given by Eq. (4.7). Besides the 
GRSVOl parton density sets we also include the sets presented in [41]. 

We show the LO double longitudinal spin asymmetry versus Pt in Fig. 
9a and the corresponding NLO result in Fig. 9b for Q = 6 GeV/c and 
= 200 GeV. They are integrated over y and the results are given in 
percent. The BB set 1 shows the most dramatic rise in All when px increases, 
while the other scenarios only show moderate increases. This is because the 
polarized gluon density in this set is larger than those contained in the other 
sets. Note that in the GRSVOl valence scenario All is negative in LO and 
NLO below pt ^ 5 GeV/c which can be inferred from Figs. 2a and 2b. 
However if the value of All(pt) is as small as indicated here this effect will 
not be measurable. A more direct comparison of the NLO results for All{pt) 
in Fig. 9b with the LO results in Fig. 9a is provided by the ratios of the 
two plots which is given in Fig. 9c. The strange behaviour of the GRSVOl 
valence scenario below = 6 GeV/c is due to the change in signs in the 
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LO and NLO results. We see that this ratio is close to unity for the GRSV 
standard scenario and the BB, set 1 while it is smaller for the BB, set 2. 
This reflects the fact that the X-factors in the polarized and unpolarized 
Pt distributions are almost equal. In Figs. 10a and 10b we show the LO 
and NLO results for All(PT)I/) as a function of y at Q = 6 GcV/c and 
Pt = 6 GcV/c. For these values the qg contribution dominates and the cross 
section is positive. Since the cross sections for the GRSVOl valence scenario 
in Figs. 4a,4b are small at this pt and fluctuate in sign across the rapidity 
range it is likely that Ai^i^{pT, y) will not be large as a function of y so we have 
omitted this parametrization in Figs. lOa-c. Again we see that the BB set 1 
distributions yield the largest results. We also take the ratios of the plots in 
Figs. 10b divided by those in Figs. 10a, which are presented in Fig. 10c. The 
latter reveals that the ratio between the LO and NLO corrected longitudinal 
asymmetry deviates considerably from unity for the BB sets in particular 
for set 2. This means that the 7^-factors for these sets differ from those 
computed for the unpolarized cross sections. In the central rapidity region 
the results are below unity reflecting the negative NLO contributions found 
earher. From Figs. 9b, 10b we can conclude that one can only distinguish 
between the various parton densities when the polarized gluon density is 
very large. In the case of the GRSV-set one cannot measure the difference 
between the valence and the standard scenario. Even the difference between 
the GRSV scenarios and scenario 2 of the BB-set is very small. If we assume 
that the proton beams have 100 % polarization (75% is more hkely) we need 
to know the polarized cross sections up to 7% to obtain Aj^l in Fig- 9b 
at Pt = 20 GeV/c up to 12.5%. This is necessary to distinguish between 
the GRSV set and BB, set 2 which is very unlikely. However if the error 
on the polarized cross sections is 17% then All at pt = 20 GeV/c can be 
determined up to 25% which is sufficient to distinguish between scenarios 1 
and 2 of the BB-set. 

Summarizing the above we have calculated the NLO corrections to the 
single particle inclusive distributions for lepton-pair production in polarized 
proton-proton scattering. We used the approximation that the lepton pair 
emerges from a virtual photon only and we neglected the additional contri- 
butions coming from the Z-boson production which is quite valid as long as 
Q < 50 GeV/c. A part of the calculation was performed using two differ- 
ent regularization methods to deal with the spurious terms introduced by 
the prescription for the 75-matrix if the matrix elements are computed in 
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n dimensions. The coefficient functions are presented in the standard MS 
scheme. In this scheme the non-singlet polarized and unpolarized coefficient 
functions satisfy the relation AWqq = —Wgq which are strong checks on our 
calculations. We have presented some NLO results for quantities of experi- 
mental interest. The fact that the qg channel yields the largest contribution 
at moderate Pt indicates that experimental measurements of the DY px dis- 
tribution will give us information on the polarized gluon density provided the 
latter is sufficiently large. We found that the K factors are moderate at least 
in the small pT-^angc and that there is a reduction in the scale dependence of 
the inclusive distributions in NLO as compared with LO. It is now up to the 
experimental groups working on polarized proton-proton scattering at RHIC 
to measure these DY distributions. 

Acknowledgement: We thank M. Klasen for discussions regarding the 
plots in [10] and C. Coriano for a clarification of the results in [14]. 
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Appendix A 

Here wc list the soft-plus-virtual gluon contributions to the polarized DY 

coefficient functions which arc determined by the one-loop corrections to the 
processes in Eqs. (2.25), (2.26) and the soft gluon corrections due to reactions 
(3.7), (3.9). The qq reaction yields the result 
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which, apart from an overall minus sign, is equal to the unpolarized soft- 
plus-virtual gluon contributions to the unpolarized DY coefficient function. 
The qg subprocess provides us with the result 
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Figure Captions 

Fig. la. The LO differential cross section dPlS.a/dQ dpr at = 200 GeV 

for the GRSVOl standard scenario plotted in the range 2 < px < 
30 GeV/c with Q = 6 GeV/c^ and /i^ = + p^^. The LO plots are 
for the subprocesses qq (dot-dashed line), qg (dashed line) and the sum 
(dotted hne). For convenience we also add the NLO result (solid hue) 
from Fig. lb. 

Fig. lb. Same as Fig. la but now for NLO. Further we have shown the con- 
tributions from all subprocesses qq (dot-dashed line), qg (long-dashed 
line), gg (dotted hne) and qq (short-dashed hne). The latter two reac- 
tions yield negative contributions we plot their absolute values. 

Fig. 2a. Same as in Fig. la but now using GRSVOl valence scenario. Since 
the qq contribution is negative at small pt we plot its absolute value. 
For convenience we also add the NLO result (sohd line) from Fig. 2b. 

Fig. 2b. Same as in Fig. lb but now using GRSVOl valence scenario. Again 
we plot absolute values when the contributions are negative. 

Fig. 3a. The LO differential cross section d^Aa/dQ dpr dy at ^/s — 200 
GeV for the GRSVOl standard scenario with Q — 6 GeV/c^, pr — 
2 GcV/c and /i^ = + Pt- The LO plots are for the subprocesses 
qq (dot-dashed line), qg (dashed hne) and the sum (dotted line). For 
convenience we also add the NLO result (solid line) from Fig. 3b. 

Fig. 3b. Same as Fig. 3a but now for NLO. Further we have shown the con- 
tributions from all subprocesses qq (dot-dashed line), qg (long-dashed 
line), gg (dotted line) and qq (short-dashed line). 

Fig. 4a. Same as in Fig. 3a but now using GRSVOl valence scenario and 
at pr = 6 GeV/c. 

Fig. 4b. Same as in Fig. 3b but now using GRSVOl valence scenario and 
at Pt = 6 GeV/c. 

Fig. 5. The dependence of d^Aa/dQ dpr on px at = 200 GeV with the 
GRSVOl standard scenario and /j,'^ — Q"^ + p^. Reading from top to 
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bottom the solid lines are the NLO results for Q = 2, 6, and 10 GeV/c^ 
and the dashed lines are the LO results. 

Fig. 6. The dependence of d^Acr/dQ dpT dy on y at = 200 GeV and 
Pt = 2 GeV/c with the GRSVOl standard scenario and /^^ = +Pt- 
Reading from top to bottom the sohd hnes are the NLO results for 
Q — 2, 6, and 10 GeV/c^ and the dashed lines are the LO results. 

Fig. 7a. The quantity N{pt,h/hq) (Eq. (4.8)) plotted in the range 0.2 < 
jJi/jJiQ < 5 (logarithmic scale) with Q — % GeV/c^ and //q = Q"^ -\- p^- 
The results are shown for px = 2 GeV/c (solid line), px = 10 GeV/c 
(dashed line), px = 20 GeV/c (dot-dashed line). The three upper 
curves on the left hand side are the LO results whereas the three lower 
curves are the NLO results. 

Fig. 7b. The quantity N{pT,y, f^/ fJ^o) (analogous to Eq. (4.8)) plotted at 
Pt — 2 GeV/c in the range 0.2 < fi/no < 5 (logarithmic scale) with 
Q — 8 GeV/c^ and nl = Q'^ + p^- The results are shown for y = ±1 
(solid line) and y — (dashed line). The two upper curves on the left 
hand side are the LO results whereas the two lower curves are the NLO 
results. 

Fig. 8a. The X-factor (Eq. (4.9)) for d'^Aa/dQ dpr plotted in the range 
2 < < 30 GeV/c for /i^ = Q"^ + p'^. The values for Q are Q = 
2 GeV/c^ (sohd line), Q = 6 GeV/c^ (dashed line) and Q = 10 GeV/c^ 
(dot-dashed line). 

Fig. 8b. The iT-factor (Eq. (4.9)) for d^Aa/dQ dpr dy at pr = 2 GcV/c 
with /j^ = + Pt- The values for Q axe Q = 2 GeV/c^ (sohd line), 
g = 6 GeV/c2 (dashed hue) and Q = 10 GeV/c^ (dot-dashed line). 

Fig. 9a. The LO longitudinal asymmetry All(j5t) (Eq. (4.7)) in % plotted 
in the range 2 < < 30 GeV/c for Q = 6 GeV/c^ and ii^ = Q"^ + 
p\. GRSVOl, standard scenario (solid hne), GRSVOl, valence scenario 
(dashed line), BB, set 1 (dot-dashed line) and BB, set 2 (dotted hne). 

Fig. 9b. Same as in Fig. 9a but now for All(pt) (Eq. (4.7)) in NLO. 

Fig. 9c. The ratio of ^ll(pt) in NLO in Fig. 9b divided by ^ll(pt) in LO 
in Fig. 9a. 
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Fig. 10a. The LO longitudinal asymmetry Ai^i^{pT,y) (Eq. (4.7)) in % at 
Pt = 6 GeV/c for g = 6 GeV/c^ and /x^ = Q^+p^. GRSVOl, standard 
scenario (solid line), BB, set 1 (dot-dashed hne) and BB, set 2 (dotted 
line) . 

Fig. 10b. Same as in Fig. 10a but now for ^LL(PT,y) (Eq. (4.7)) in NLO. 

Fig. 10c. The ratio of All(pt, y) in NLO in Fig. 10b divided by All(pt, v) 
in LO in Fig. 10a. 
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